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Quantum Noises, Physical Realizability and Coherent Quantum 

Feedback Control 

Shanon L. Vuglar, Ian R. Petersen 


Abstract —Physical Realizability addresses the question of 
whether it is possible to implement a given linear time invariant 
(LTI) system as a quantum system. A given synthesized qnantnm 
controller described by a set of stochastic differential equations 
does not necessarily correspond to a physically meaningful quan¬ 
tum system. However, if additional quantum noises are permitted 
in the implementation, it is always possible to implement an 
arbitrary LTI system as a quantum system. In this paper, we 
give an expression for the number of introduced noise channels 
reqnired to implement a given LTI system as a quantum system. 
We then consider the special case where only the transfer function 
to be Implemented is of Interest. We give results showing when it 
is possible to implement a transfer function as a quantum system 
by introducing the same number of quantum noises as there are 
system outputs. Finally, we demonstrate the utility of these results 
by providing an algorithm for obtaining a suboptimal solution 
to a coherent quantum LQG control problem. 


I. Introduction 

For systems where it is necessary to consider quantum 
effects, the laws of quantum mechanics introduce new consid¬ 
erations not present in classical controller synthesis problems. 
The presence of quantum noises ID introduce fundamental 
limits on controller performance. Furthermore, the requirement 
for unitary evolution, the non-commuting nature of quantum 
observables, and the requirement for commutation relations to 
be preserved as systems evolve (see for example ID) lead to 
the notion of physical realizability El, 0. This is the property 
that a given system model represents the dynamics of a phys¬ 
ically meaningful quantum system. Controller synthesis and 
optimization problems that are well understood in the classical 
regime can become difficult when restricting their solutions to 
physically realizable quantum controllers. New and tractable 
methods are required for these quantum controller synthesis 
problems. 

It is useful to draw a distinction between measurement based 
quantum feedback control and coherent quantum feedback 
control. In measurement based quantum feedback control, 
measurements of observables of a quantum system are used 
to apply feedback via a classical controller. While the closed 
loop system is modeled and analyzed in a quantum setting, 
the results of the measurements are classical signals and the 
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controller can be implemented using analog or digital elec¬ 
tronics. This paper addresses the alternative coherent quantum 
feedback approach illustrated in Figure [T] in which quantum 
systems are interconnected directly, avoiding measurement; 
e.g., see 13. 
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Fig. 1. A Coherent Quantum Control Scheme 


Both forms of quantum feedback control are relevant to a 
diverse range of applications which take advantage of quan¬ 
tum effects. These applications include quantum computing, 
quantum communications, quantum cryptography and preci¬ 
sion metrology such as gravity wave detection. Measurement 
based quantum feedback control is well understood (see for 
example 0) and has been used successfully to manipulate 
quantum effects. For example, in Q a qubit was maintained 
in an oscillating superposition state. This important result is 
relevant to the field of quantum computing. 

The majority of experimental results to date have focused 
on measurement based feedback control. Coherent quantum 
feedback control presents additional challenges in that the 
controller must be physically realizable. However, coherent 
quantum feedback control may offer several advantages over 
its measurement based counterpart. Firstly, coherent feedback 
avoids the collapse of the quantum state and the loss of 
quantum information associated with the use of quantum mea¬ 
surement. This is particularly relevant to quantum computing 
where quantum states need to be maintained and manipulated. 
Secondly, implementing coherent controllers may introduce 
fewer quantum noise channels than the measurement process 
and this in turn may lead to better control system performance. 
Finally, it may be that there are technical or experimental 
benefits in implementing a controller as a quantum system. 
For example, the use of a coherent controller may result 
in advantages in terms of the speed of control. Also, the 
experimental setup may make measurement impractical. 

Coherent quantum feedback control is generating increasing 
interest within the research community ( El, El, El, El) and 
central to this area is the notion of physical realizability. In 
the classical setting, we regard controllers as always being 
possible to implement. In the quantum setting, a given syn¬ 
thesized quantum controller described by a set of stochastic 
differential equations cannot always be implemented by a 
physically meaningful quantum system. Several recent papers 





have addressed this issue of physical realizability 0,11,111, 
d, giving conditions for when a given system is physically 
realizable. Other papers ifT^ - lfTsl have given algorithms for 
experimentally implementing several classes of physically 
realizable quantum systems. 

An importance difference between classical and quantum 
controller synthesis is that in the case of coherent quantum 
feedback control, implementing a controller as a quantum 
system may require the introduction of quantum vacuum 
noises. To see how this might arise, consider the following 
example from quantum optics. 
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Fig. 2. An Optical Cavity 

Suppose that as part of the controller implementation pro¬ 
cess, the design calls for a laser beam to be passed through 
an optical cavity as shown in Figured Here a naive approach 
would be to consider this device as having a single input and 
single output. However, the laws of quantum mechanics imply 
that there is a second input to the cavity. Indeed, the mirror on 
the right, which produces the output, also causes the cavity to 
be coupled to a vacuum noise input; e.g. IH. To obtain correct 
results when modeling such a system, it is essential to take 
this additional vacuum noise source into account. 

Utilizing well established controller synthesis methods 
(such as H°° controller synthesis) and modifying the classical 
solutions by incorporating additional quantum noises to obtain 
physically realizable quantum systems provides a tractable 
approach to coherent quantum controller design; e.g., see 0. 
This approach requires a method for determining how many 
additional quantum noises are necessary for physical realiz¬ 
ability, and for constructing the resulting quantum systems. 

In 0, the authors demonstrated that it is always pos¬ 
sible to implement an arbitrary, strictly proper, linear time 
invariant (LTI) system as a quantum system by introducing 
a sufficient number of quantum vacuum noise channels. It is 
straightforward to obtain upper and lower bounds on the num¬ 
ber of introduced quantum noises that are necessary to obtain 
physical realizability. Since these noises place limits on the 
achievable controller performance, it is desirable to minimize 
the number of these introduced noises. This paper extends the 
result in 0 to determine the number of introduced quantum 
noises that are necessary to implement a given, strictly proper, 
LTI system. Also, our result extends the construction method 
in 0 to give a construction that only introduces as many 
quantum vacuum noises as are necessary to make that system 
physically realizable. 

We also consider the special case in which we are only in¬ 
terested in physically realizing a transfer function, as opposed 
to a specific state space realization. Note that the number of 
introduced quantum noises necessary to physically realize a 


strictly proper LTI system must be at least as many as the 
output dimension. We provide a condition, under which a 
strictly proper transfer function can be physically realized with 
the number of introduced quantum noises being equal to the 
output dimension. This condition is given in terms of a non¬ 
standard algebraic Riccati equation. We then provide condi¬ 
tions for the existence of a suitable solution to this Riccati 
equation. This leads to a numerical solution to the question 
of whether a particular strictly proper transfer function is 
physically realizable with the number of introduced quantum 
noises being equal to the output dimension. 

Preliminary conference versions of the results of this paper 
have appeared in ifThll - lfTSl . Here, we provide detailed proofs 
not included in those conference papers. Furthermore, we 
demonstrate the utility of our main results by providing an 
algorithm to obtain a suboptimal solution to a coherent quan¬ 
tum linear quadratic Gaussian (LQG) problem. This algorithm, 
and the example demonstrating its application, did not appear 
in the conference papers. 

The remainder of the paper proceeds as follows. In Sec¬ 
tion m we describe the quantum systems considered through¬ 
out this paper. In Section [Till we recall the definition of 
physically realizable systems and outline relevant previous 
results on this topic. In Section |IV] we present our main 
results. We first consider the problem of implementing a 
particular state space model as a quantum system and the 
number of introduced quantum noises necessary to do so. We 
then consider the special case where a transfer function is to 
be physically realized. We give results regarding when such a 
transfer function is physically realizable with the number of 
introduced noises being equal to the dimension of the system 
output. In Section |V] we demonstrate the utility of our results 
by presenting an algorithm for finding a quantum controller 
which is a suboptimal solution to a coherent quantum LQG 
problem. An example demonstrating our algorithm followed 
by our conclusion are then given in Sections |VT] and IVIII 
respectively. 


H. Quantum Systems 
A. General Quantum System Model 

Open quantum harmonic oscillators represent an important 
class of quantum systems. Such systems can be described 
by quantum stochastic differential equations (QSDEs) of the 
following form (see 0): 

dx{t) = Ax{t) dt + B dw(t); 

dy(t) = Cx{t) dt + D dw(t). (1) 

Here, X{t) = [xi{t) ••• Xnit)]'^ is a column vector 
of n self-adjoint system variables which are operators on an 
underlying Hilbert space. Being quantum in nature, these vari¬ 
ables do not commute in general. The commutation relations 
for these variables are described by a real skew-symmetric 
matrix 0: 

[xi{t),Xj{t)] = Xi{t)xj{t) - Xj{t)xi{t) = 2iQij. 

Similarly, dw{t) is a column vector of riw self-adjoint, non- 
commutative operators representing the input to the system and 









dy{t) is a column vector of Uy self-adjoint, non-commutative 
operators representing the output of the system. Their com¬ 
mutation relations are as follows: 

[dtyi(f), dwj(f)] = 

[dy2(f), dyj(f)j — 2z0y df, 

where and Qy are real skew symmetric matrices. 

The input signals dw{t) are assumed to admit the decom¬ 
position 

dw{t) = /?u,(f) df -I- dui(f) 

where the self-adjoint, adapted process /3^(f) is the signal 
part of dw{t) and dw{t) is the noise part of dw{t). Here, 
Pw{t) is assumed to commute with x{t). The vector dw{t) is 
a quantum Wiener process with Ito products 

dw{t) dw'^{t) = dt 

where is a non-negative Hermitian matrix. Let F^ = 
Sw + TiD, where is real and is imaginary. Then 
describes the intensity of the quantum Wiener process 
and is the quantum analog of the intensity matrix for a 
classical Wiener process. The commutation relations for dw 
are determined by T^: 

[dti;(f), dw^jf)] = 2Twdt. 

Since /3u,(f) is an adapted process, /?u,(f) commutes with 
dw^t) for all t > 0. Also, dw{t) commutes with x{t). 

Finally, n, and Uy are even (this is because in the quan¬ 
tum harmonic oscillator, the system variables always occur as 
conjugate pairs, see 0) and A, B, C and D are appropriately 
dimensioned real matrices describing the dynamics of the 
system. For further details regarding these models, see El. 

Remark 1: While it is always possible to describe a col¬ 
lection of quantum harmonic oscillators by QSDEs of the 
form O- not all QSDEs of this form correspond to a col¬ 
lection of quantum harmonic oscillators. The property that the 
QSDEs ([Til correspond to a collection of quantum harmonic 
oscillators is called physical realizability and is addressed in 
greater detail in Section [III] 

We will now further restrict our attention within the class 
of quantum systems described above. 


B. A Class of Quantum System Models 

This paper addresses the problem of implementing an 
arbitrary, strictly proper, LTI system as a quantum system 
(for example when implementing a coherent controller) by 
introducing vacuum noise sources. The resulting quantum 
systems are described by the following QSDEs which are a 
special case of dD: 

dx{f) = Ax{t) dt + Bu du{t) 

+ By^ dvi{t) + By^ dv 2 it); 
dy{t) = Cx{t) dt + dvi{t). (2) 

Here, du{t) (a column vector with Uy components) represents 
the inputs to the system and, like d'w{t) in ([T]), is assumed 
to admit the decomposition dM(f) = f]u{t)dt+ du{t). Also, 


dvi{t) and dv 2 {t) (column vectors with and Uy^ compo¬ 
nents respectively) are quantum Wiener processes correspond¬ 
ing to the introduced vacuum noise inputs. Eor convenience, 
the vacuum noises are partitioned into two vectors dui (f) and 
dv 2 {t) such that = n„. Then, Uy = Uy^ + Uy^ is the total 
number of introduced vacuum noise inputs. Subsequently, we 
will refer to dui as the direct feedthrough quantum noises 
and to dt ;2 as the additional quantum noises. Also, Fy, Sy, 
Ty, Fy^, Sy^, Ty^, Fy^, Sy^ , and Ty^ aK defined for du{t), 
dui(f) and du 2 (f) respectively as Fjj, and were for 
d?i;(f) in ([T]). Eurthermore, we assume that Fyi and Fy 2 are 
appropriately dimensioned block diagonal matrices with each 
diagonal block equal to [ j ]. This assumption corresponds 
to the fact that dui and dv2 represent vacuum noises El- 
The remaining symbols have the same meanings as in ([T]). We 
restrict our attention to the case where Uy = riy. 


HI. Physical Realizability 

A. Definitions 

As in El, 0 , cni, M-M, the concept of physical 
realizability means that the system dynamics described by 
the QSDEs ([T]) correspond to those of a collection of open 
quantum harmonic oscillators. Here, we slightly modify the 
definition of physically realizable given in El Definition 3.1]. 
In 0 both fully quantum systems and hybrid systems with 
quantum and classical degrees of freedom are considered. 
However, we restrict our definition of physically realizable 
systems to those that are fully quantum. 

Definition 1: The system variables x are said to satisfy the 
canonical commutation relations if 




[ 

Xi{t) 

Xj{t)\ = 2iQij 




where 0 is of the form 
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of pairs of position and momentum operators: 
[qi{t) pi{t) 92(f) P2{t) ...]^. 

Definition 2: The system described by dB is physically 
realizable if 0 is of the form 0 and there exists a quadratic 
Hamiltonian operator B = |x(0)^i?x(0), where i? is a real, 
symmetric, n x n matrix, and a coupling operator vector 
£ = Aa;(0), where A is a complex-valued iriu, x n coupling 
matrix such that the matrices A, B, C and D are given by: 


A = 2Q{R + '3m (AU)) ; 
B = 2iQ [-At A"^] F; 


Sn„ 0 

A + A* 


—iA + iA* 


D= [In 


0 


Uy X {n, 


,-ny)\ ■ 


(4a) 

(4b) 

(4c) 

(4d) 


Here: = Pdiag(M); M = \[\Ii\\ 

^ny = l.^k'^yX^ny ; P is 










the appropriately dimensioned square permuta¬ 
tion matrix such that P [ai 02 • • • a 2 m] 

= [ai as ■■■ a 2 m-ia 2 04 ■ ■ ■ a 2 m] and diag(M) 
is an appropriately dimensioned square block diagonal matrix 
with each diagonal block equal to the matrix M. (Note that 
the dimensions of P and diag(M) can always be determined 
from the context in which they appear.) 3m(.) denotes 
the imaginary part of a matrix and ^ denotes the complex 
conjugate transpose of a matrix. 

Remark 2: This definition amounts to saying that a sys¬ 
tem ([T]i is physically realizable if and only if it corresponds 
to a collection of open quantum harmonic oscillators 13. 

We now apply this definition to the class of quantum 
systems (H. The system (|2]l is physically realizable if there 
exists a real, symmetric, nxn matrix R, and a complex-valued 
i + n „2 -f n„) X n matrix A such that the matrices A, 

By^, By^ and C are given by: 


[B 


V\ Pv2 


A = 2Q{R + 3m (AU)) ; 
By] = 2i& [-At A^; 


C = P" 


'Briy 0 


A -1- A?^ 

0 


-ik + iA* 


(5a) 

(5b) 

(5c) 


where 0 is of the form Q- Here, Y,ny = 

^riyX^riy ^inyX^^ny^+ny^+Tiy-riy') ■ 

Theorem 1: (see 13 Theorem 3.4]) The system ([T]i is phys¬ 
ically realizable if and only if: 


0 = iAQ -f iQA^ + BT^B ^; 


B 


= 0C^diag(J); 


where 0 is defined as in ([3 and D satisfies (l4dl i. 

Corollary 1: The system (|3 is physically realizable if and 
only if: 


0 = iAQ + iQA^ + By, Ty, Bi, 


TD rp tdI 

^V2-^V2-^V2 


By 


= 0C^diag(J), 


where 0 is defined as in (13- 


BuTuBy ; 


B. Previous Results 

In 13, it was demonstrated that by introducing a sufficient 
number of vacuum noises, an arbitrary LTI system could be 
made physically realizable. In particular, the following lemma 
relating to the physical realizability of a purely quantum 
controller was proved. 

Lemma 1: (See (3 Lemma 5.6]). Let Fy be a block di¬ 
agonal matrix with each diagonal block equal to I + iJ, 
and let A,B and C be matrices such that A € S 

]^nxn„ g jg( 0 defined as in ([ 3 - 

Then there exists an even integer Uy^ > 0 and matrices 
By, G g ]^nxra„ 2 ^ system (13 is 

physically realizable. 

Remark 3: It follows from 13 Theorem 3.4], that for a 
system described by a strictly proper transfer function 

G{s) = C {si - A)-^ By, 


the dimension of the system output Uy is a lower bound on 
the total number of introduced vacuum noises Uy that are nec¬ 
essary for the system to be physically realizable. That is, the 
direct feedthrough quantum noises dui in the system ([3 are 
necessary, but may not be sufficient for physical realizability. 
We are also interested in the situation in which the presence 
of the noises dui is sufficient for physical realizability and 
the noises dt ;2 are not needed. In this case, we say that the 
LTI system is physically realizable with no additional vacuum 
noises. Physically realizing a system with minimal additional 
noises means to implement the system as a quantum system (|3 
by only introducing the number of additional vacuum noises 
Uy^ > 0 , that are necessary for physical realizability. 


IV. Main Result 

A. General Case - Implementing a State Space representation 

In this section, we give a method to physically realize a 
strictly proper LTI system 

da; = Ax df -I- By du 
dy = Cx dt 


with minimal additional quantum noises. The remainder of this 
section is structured as follows. We first give our algorithm. 
We then formally state our result as a theorem. The subsequent 
proof of the theorem justifies our algorithm. 

The algorithm for obtaining a physically realizable sys¬ 
tem (13 with minimal additional quantum noises proceeds as 
follows: 

1) Construct the matrix 

5 = 0P„0„„Pj0-0A-A^0-C^0„„C. (6) 


Here 0„^ and 0„ are commutation matrices of the 
form (13 of dimensions Uy x n„ and Uy x riy respectively. 


2 ) 


3) 

4) 

5) 

6 ) 


Find the rank of the matrix S: r 


rank 


S' 


. Now 


ny„ = r. That is, direct feedthrough quantum noises, 
and r additional quantum noises, are necessary for the ex¬ 
istence of By, , By^ such that the system (13 is physically 
realizable. This gives Uy = Uy + r. 

Calculate S = |S. 

Construct the singular value decomposition (SVD) for S: 
S = DU. Here D is diagonal and U is unitary. 
Construct A{,i = (|D| + D)^ U where \D\ is the diagonal 
matrix with entries equal to the absolute values of the 
corresponding entries in D. 

Construct By, and By^ as follows: 


By, = 0C'^diag(J); 

By^ = 2i0 [-A|[^ A^JPdiag(M). 

The system (13 with {A,By,C} given and with By,, By^ so 
constructed is physically realizable. We now give a theorem 
which formally states that Uy = Uy + r introduced noises 
are necessary and sufficient for physical realizability. The 
construction of By, and By^ above follows from the proof 
of the theorem. 

Theorem 2: Consider a strictly proper LTI system defined 
by given matrices A, By and C. There exist matrices By, 












and such that the corresponding system ([2]) is physically 
realizable and with equal to r where r is the rank of 
the matrix 0 — 0A — A^0 — C^0„j^C'). Con¬ 

versely, suppose that there exist matrices and such 
that the corresponding system (l2]i is physically realizable. Then 
n „2 > r. 

Proof: The proof is structured as follows. We first show 
that n„ + r introduced quantum noises are sufficient for 
physical realizability. We then show that + r introduced 
quantum noises are necessary for physical realizability. 

Following the method of lO, the construction of the matri¬ 
ces R, A, By^ and By^ in (fSal) - (l5^ is as follows: 


R 

A 


By^ 

By2 


-i(0A + A^0^); 




2 



0C^ diag(J); 



2*0 [-Aj, A^JPdiag(M), 


where 0 is defined as in (O. Here, 


Ab 2 = -i[l 0] PdisLg{M)ByQ; 
and Abi is any complex ^ny^ x n matrix such that 

Aj^iAfoi = L^i 


/A^0^ — 0A IpTpT 

■ 0 1 

V 4 4^ ^ 

0 

1 

_ 1 


— 3m 


(M>2^b2 


PC 

(7) 


where Si is any real symmetric nxn matrix such that A|^j^A^ i 
is nonnegative definite. 

The matrix At,i can be constructed as follows: first a real 
symmetric nxn matrix Si is constructed such that the 
right hand side of ([Til is nonnegative definite. Then A{,i is 
constructed such that (|7]i holds. 

Note that A^i has ^ny^ rows and thus determines the 
number of additional quantum noises required in this imple¬ 
mentation. We now provide a method for choosing Si and 
Afii to obtain the result. 

It is desired to construct Si such that A^^^^Af,! is of minimum 
rank. This will allow A{,i to be constructed with the minimum 
number of rows. We make the following observations about 
the terms in equation (|7]i: 




- 0 
0 /] p 

-I oj 

^h2^b2 


3m (a]; 2 A 62 ) 


= e^; 

= 0n„; 

= 0^p.rt 


( 8 ) 

(9) 


[/ ojrs;:© 


--0P„ diag 
-^QBu diag 


1 i 
—i 1 
0 1 
-1 0 




B.Q 


= --0P„0„„P^0. 


( 10 ) 


Substituting ([8]l, (|9]) and (fTOl i into (|7]i, we obtain 


where 

^2 = 

and S is defined as in (|6]l. 

Note that the matrix S is real and skew symmetric. Thus 
S = m Hermitian, has real eigenvalues and is diagonaliz- 
able: S = U^DU where D is diagonal and U is unitary. 

We wish to find a real, symmetric matrix Si such that S 2 = 
Si -f S' is positive semi-definite and of minimum rank. Let 

51 = C/1 |i3| t/. We claim that Si is real and symmetric, that 

5 2 = Si -f S > 0, and that S 2 has rank equal to half that of 
S. 

First, we show that this matrix Si is real and symmetric. 
Observe that Si = Sil and Si > 0. Also: 

Si^ = c/i|p|V = c/tp^c/ = S^ 

Here, S is purely imaginary, thus S^ is real and Si^ > 0 is also 
real, and therefore has a real square root. From the uniqueness 
of the positive semi-definite square root of a positive semi- 
definite matrix na Theorem 7.2.6] we conclude that .^1 is 
real. Further, since Si is Hermitian, Si is symmetric. 

We now show that S 2 has rank equal to half that of S 
and that S 2 is positive semi-definite. We observe that S is 
Hermitian and so its eigenvalues are real. Thus the eigenvalues 
of S are purely imaginary. Also since S is real, its eigenvalues 
occur in complex conjugate pairs. Thus D is of the form: 


Ai 0 0 0 

0 -Ai 0 0 

0 0 A 2 0 

0 0 0 -A 2 


; A, > 0. 


\D\ 


Ai 0 0 0 
0 Ai 0 0 
0 0 A 2 0 
0 0 0 A 2 


;A, > 0. 


\D\+D 


'2Ai 0 0 0 

0 0 0 0 

0 0 2 A 2 0 

0 0 0 0 


; A, > 0. 


From this, it can be seen that \D\+ D has a rank which is 
half that of D. Since 

^2 = ^1 + 5" 

= C/+|i3|[/-f C/^PC/ 

= uWd\ + d)u, 

it follows that S 2 is positive semi-definite and has a rank which 
is half that of S. 
























Since S and S have the same rank, S 2 has rank | where 
r the rank of S. Since S 2 > 0 has rank it is possible to 
construct A^i with | rows, such that S 2 = Recall 

that, Af,i has rows, and we have riy^ = r. That is, the 
system is physically realizable with the number of additional 
quantum noises riy equal to + r where r is the rank of the 
matrix S defined in (|6]l. 

We now consider the second part of the theorem and show 
that + r introduced noises are necessary for physical 
realizability. To do so, it is sufficient to show that the number 
of columns of must be greater than or equal to 

n„ + r. 

From (fSbI i. it can be shown that: 


Bu 

= 2iQ 

l-A 

^ 12 ] r; 

(lla) 

By^ 

= 2iQ 

[—^60 

A^o] r; 

(lib) 

By^ 

= 2i0 

[-Ail 

AfJ r; 

(11c) 


where 


A = 


Af,o 

Afoi 

Ab2 


That is, By^ has twice the number of columns as A^i has rows. 
Therefore, we wish to show that A^i has at least ^ rows. 
Consider, 


3m{A^A) = 3m(A^pAbo) + Tm(A|,^Abi) + 3m{Al^Ab2). 


That is, 

3m(A^;^At,i) = 3m(A'fA)-3m(AjgAf,o)-3m(Aj2A62)- (12) 

Rearranging (l4al i. we obtain 

= R + 3Tn(AlA), 

where R and 3m(AlA), are respectively the symmetric and 
skew-symmetric parts of the left hand side of this equation. 
From this, it can be shown that 

3m(AlA) = -^(0A -f A^0). (13) 

Also using (l4cl i and (|9]l, it is straightforward to verify that 

Tm(AJoAbo) = Ic^QuyC. (14) 

Substituting (fTOl i. (fTSl l and (fT4l) into (fT^ we obtain 
Tm(AjiAbi) = 


where S is defined as in That is. 


A^liAbi — -b 


where Si is the real part of A^^Af,i. 


Now using II 20 I Fact 2.17.3], we observe that 


rank 


(aLa.i) 


= rank 


= -rank 
2 


> -rank 
- 2 


> -rank 
- 2 


That is, for any 


rank 


{M 


blAbl 


+ I — 


" S 
4 
S 

.“4 “1. 

" S 
-1 4 


> -rank]^]. 


This in turn implies that Abi has at least ^ rows, where, r 
is the rank of the matrix S defined as in (ISl). However, it 
follows from (fTlB . that has twice as many columns as 
Abi has rows. That is, By^ has at least r columns and hence 
[By^ 5 ^ 2 ] has at least Uy + r columns. Hence, the number 
of quantum noises ny is greater than or equal to + r. This 
concludes the proof of the theorem. ■ 


B. Special Case - Physically Realizing a Transfer Function 
When designing LTI controllers, usually the transfer func¬ 
tion of the controller rather than its particular state space real¬ 
ization determines the closed loop performance. As such, the 
question of whether a particular transfer function is physically 
realizable may be of greater interest than whether a particular 
state space realization is physically realizable. 

Therefore, we now turn our attention to the case in which 
we are interested in implementing an LTI quantum system 
with a specified strictly proper transfer function. This is 
equivalent to allowing a state transformation on the state 
space description of the system. In particular, we consider 
the problem of whether a particular transfer function can be 
physically realized by only introducing the direct feedthrough 
quantum noises dui in (|2]). That is, without introducing any 
additional quantum noises dv 2 - 

Here, we recall from Remark [3 that for systems described 
by strictly proper transfer functions the direct feedthrough 
quantum noises are necessary for physical realizability. 

Under some assumptions, it is possible to implement a 
specified transfer function as a physically realizable quantum 
system (|2l) where only the direct feedthrough quantum noises 
are introduced. 

Theorem 3: Consider a system with strictly proper transfer 
function matrix: 

G{s) = C{sI-A)-^Bu. 

Suppose the algebraic Riccati equation (ARE) 

XByQy^B^X -A'^X-XA-C^en^C = 0 (15) 

has a non-singular, real, skew-symmetric solution X. Here, 
the matrices 0„^ and 0„^ are defined as in (|3ll. Then there 
exists matrices {A,Bu,C} such that 


G{s) = C{sI-A)-^By 










and the corresponding system 0 is physically realizable 
with only the direct feedthrough quantum noises dui and no 
additional quantum noises dt; 2 . 

Proof: First, note that for any 2m x 2m non-singular 
real, skew-symmetric matrix X there exists a non-singular, 
real matrix T for which X = T^0T where 0 is defined as 
in © lEU Corollary 8.24], Let, 

X = T G R”; det T ^ 0; 

A = TAT-^; 

Bu = TBu\ 

C = CT-i; 

= 0C^diag(J). 

The result now follows by applying Theorem|2] Indeed, (flSl l 
implies; 

0 = XBuQn^BlX-X^X-XA-C^Qn^C 

= T^QTlBu&n^BlT'^QT - A^T^QT-T'^QTA 

-C^On^C 

= QTBuQu^bIt'^Q - (T^)”^ i^r^0 - STAT-^ 

- {T^y" C^Sn.CT-^ 

= 0TB„ 0„„ (tBu^ ^ 0 - (tAT-^ Y e- QTAT-^ 

- (cT-^y Qn,CT-^ 

= QBuQn^BjQ- A^Q-QA-C^Q^^C. 


That is, the matrix S defined in ® has rank zero. Applying 
Theorem m we conclude that the system {A,Bu,C} can be 
physically realized with riy^ =0. ■ 

We now give conditions for when the ARE (fTSI l has a non¬ 
singular, real, skew symmetric solution X. The proof given 
below closely follows that in E2\ Theorem 13.5]. This result 
also leads to a numerical procedure for physically realizing a 
strictly proper transfer function with the minimal number of 
additional quantum noises. 

For convenience, we define R = —BQn^B^^ Q = 
and rewrite ( fTsT i: 

A^X + XA + XRX + Q = 0. (16) 


Note that Q and R are skew symmetric. 
Define 


H = 



R 

_^T > 


(17) 


Im [y^] where Xi,X 2 G Here Im denotes the 

subspace spanned by the columns of the matrix [ Xs ] ■ 

We assume Xi is non-singular or equivalently that x-{H) 
and Im[j] are complementary subspaces. Then define X = 
X 2 Xf^. It follows that X is uniquely determined by H. 
We will denote the corresponding function by X = Ric{H) 
with the domain dom(Ric) consisting of matrices H of the 
form dnii satisfying the properties that H has no purely 
imaginary eigenvalues, and that Xi is non-singular. 

Theorem 4: Suppose H G dom(Ric) and X = Ric{H). 
Then X is skew-symmetric and solves the algebraic Riccati 
equation X + XA + XRX + Q = 0. 

Proof: Let Xi, X 2 be as above. There exists a Hurwitz 
matrix i7_ G such that 


X' 


Xi' 

to 


.^2, 


Pre-multiply (fTsT i by [ ^ Z to obtain 



Since ZH is skew-symmetric, so are both sides of the above 
equation. From the right-hand side: 

(XjXi + XfX2)iL_ = -i7l(XjXi + XfX2)^. 

This is a Lyapanov equation. Since is Hurwitz, the 
unique solution is XjXi -f X^ X 2 = 0. That is, XjX 2 
is skew symmetric, and since Xi is non-singular, X = 
{Xyy [X^ X 2 )Xf"^ is also skew-symmetric. 

It remains to be shown that X is a solution to (fThl l. Post- 
multiplying (fTsT l by Xf^, we obtain 

H[A] = [yx,H.xy 

and pre-multiplying by [x -/] gives 

{^-AH[y = 0 : 

which is precisely (fThl l. ■ 

Remark 4: The above proof also leads to a numerical 
procedure for solving the ARE (fTSi) and hence solving the 
physical realizability problem under consideration. This nu¬ 
merical procedure involves solving the eigenvalue, eigenvector 
problem for the matrix H. The following corollary, which 
follows directly from combining Theorems [3 and @1 is the 
main result of this subsection. 

Corollary 2: Consider a system with strictly proper transfer 
function matrix: 

G{s) = C{sI-A)-^Bu. 


and Z = —*[/q]. Note that Z~^ = Z\ {ZH) is skew 
symmetric, and Z~^HZ = Z^HZ = —Hf That is, H and 
—W are similar, from which it follows that A is an eigenvalue 
of H if and only if —A is. That is, the eigenvalues of H are 
symmetric about the imaginary axis. 

Assume H has no eigenvalues on the imaginary axis and 
let x-{H) be the n-dimensional spectral subspace Il22l of H 
corresponding to its negative eigenvalues. We find a set of 
basis vectors for x-{H) and stack the basis vectors to form 
a matrix. Partitioning this matrix, we can write X-{H) = 


Suppose H G dom(Ric) and X = Ric{H) is non-singular 
where H is defined as in dnii. Then there exists matrices 
{A, Bu,C} such that 

G{s) = C{sI - A)-^Bu 

and the corresponding system @ is physically realizable 
with only the direct feedthrough quantum noises dui and no 
additional quantum noises du 2 . 

We now give the algorithm for solving the ARE (fTSI) and 
hence physically realizing a given transfer function by only 








introducing direct feedthrough quantum noises. Suppose we 
wish to physically realize the transfer function 

G(s) = C{sI-A)-^Bu- 

1) Construct the matrix H as in (fTTI i. Find the eigenvalues 
and eigenvectors of H. Check that H has no purely 
imaginary eigenvalues. In practice, this means checking 
that the real part of each eigenvalue has magnitude greater 
than some small numerical tolerance. 

2) Construct a matrix [ ] such that its columns are the 

eigenvectors of H that correspond to eigenvalues with 
negative real part. Check that Xi and X 2 are non-singular 
and calculate X = X 2 Xi^. The matrix is a non¬ 
singular solution to the ARE (flsT l. 

3) Find the eigenvalues and eigenvectors of X. These will 
occur in complex conjugate pairs. Hence, construct a 
diagonal matrix A with entries being the eigenvalues of 
X and with complex conjugate eigenvalues in adjacent 
columns. Construct a matrix V with columns being the 
corresponding eigenvectors of X normalized to have 
unit norm and with complex conjugate eigenvectors in 
adjacent columns. 

4) Construct the n x n diagonal matrix A with alternating 

diagonal entries i and —i. Also construct the nxn block 
diagonal matrix V with each diagonal block correspond¬ 
ing to ^ ]. 

5) Calculate T = VDV^ where D = ^A“^A^ ^. 

6) Construct 

A 
Bu 
C 

Bvi 

Then 

G{s) = C{sI - A)-^Bu 

and the system (|2|i corresponding to | A, i?„, cj is physically 
realizable by introducing only direct feedthrough quantum 
noises with By^ constructed as above. No additional quantum 
noises are necessary for physical realizability. 

Remark 5: We now justify the above numerical algorithm 
for constructing T such that X = T^QT. 

Since X is skew-symmetric, X = EAE^, where E is a 
unitary matrix with columns which are the eigenvectors of 
X. Also, A is a diagonal matrix where the diagonal elements 
are the eigenvalues of X, which are purely imaginary and 
occur in complex conjugate pairs. For every eigenvector v 
of X corresponding to eigenvalue A, its complex conjugate 

V is also an eigenvector and corresponds to A. If necessary, 
we reorder the columns of V and corresponding entries of 
A such that these complex conjugate pairs are adjacent: 

V = [vi V2 V2 ... ] . 

Similarly we can write 0 = VAV^ where E is a block 
diagonal matrix with repeated blocks [J and A is a 
diagonal matrix with alternating entries i and —i. Observe 


that there exists a diagonal matrix 


di 

0 

0 

0 • ••' 

0 

di 

0 

0 ••• 

0 

0 

^2 

0 •• • 

■ 0 

0 

0 

d2 ■ ■ ■ 






such that A = DAD and the diagonal elements di are real 
and positive. We now have X = VDV^QVDV^ . Define T = 
EDE^ then X = T'^ST. 

Observe that the matrices [ J ] and [ q’ ] commute. 
Therefore E and D commute. We claim that the matrix 
T = DEE^ is real. This follows since 

(EE^')^ = EE^ 

= [«! Cl «2 C 2 ... ] diag ( ^ J “* ]) 

= [(-Ul-I-Cl) -i(vi-vi) (V2+V2) -i(v2-V2) ... ] 

which is real. Therefore T as constructed above is real and 
X = T^0T. 

V. A SUBOPTIMAL COHERENT QUANTUM LQG 
CONTROLLER DESIGN ALGORITHM 

In this section, we use the results in Section |IV] to provide 
an algorithm for the design of a suboptimal coherent quantum 
LQG controller. The main idea of our algorithm is to design a 
classical LQG controller and then use Theorem 2 to implement 
this controller as a quantum system. In contrast with the clas¬ 
sical LQG controller synthesis problem, here the separation 
principle of using the optimal state estimator and optimal 
regulator no longer applies due to the relation between the 
optimal regulator gain and the additional quantum noises that 
arise when implementing the controller as a quantum system. 
This is addressed in the algorithm that follows. 

A. Problem Formulation 

We will now formally state the problem to be addressed. Our 
formulation follows that in 0 with some minor differences. 

Suppose we have a quantum plant described by the follow¬ 
ing QSDEs which are a special case of O: 

dx{t) = Ax{t) dt + Budu{t) + Byj^ dwi{t)-, 
dy{t) = Cx{t) dt + Dudu{t) + Dyj^ dwi{t). (19) 

Here, as in O, X is a column vector of n self-adjoint 
system variables. The column vector du{t) represents the 
input to the system. It consists of signals of the form 
du{t) = I3u{t) dt -f du{t) where du{t) is the noise part 
of du(f) with Ito products du{t)dvA{t) = Fydt where 
Fy is non-negative Hermitian. Also, the self-adjoint, adapted 
process /3„(f) is the signal part of du{t). Furthermore, dtui 
is a column vector of riwi non-commutative quantum Wiener 
processes with Ito products dtyi(f) dwi^(f) = dt where 
Fyj^ is non-negative Hermitian. Here, drui represents noises 
driving the system and may for example include vacuum 
noises and/or thermal noises. The column vector of Uy signals 
dy{t) represents the output of the system. Finally, as in ([T]i, n. 


= 

= TBy-, 

= CT-^; 

= 0G^diag(J). 





nu, n-wi and Uy are all assumed to be even and A, B^, B^-^, 
C, Du and Du,^ are appropriately dimensioned real matrices 
describing the dynamics of the system. For further details see 
0 , 0. For simplicity we restrict our attention to the case 
where Uy = n„. 

Furthermore, suppose that we wish to minimize an infinite 
horizon quadratic cost function: 


1 

J = lim sup — / 
tf^ao tf Jq 


( 20 ) 


where (.) denotes the quantum expectation; e.g., see 0. 

We restrict attention to controllers described by the follow¬ 
ing QSDEs which are of the form (|2|i: 


dxK{t) = AKXK{t)dt + Bydy{t) 

+ Bu^ dvi{t) + By^ dv 2 (t); 
du(t) = CKXK(t) dt + dvi(t). (21) 

The problem is to design a physically realizable quantum 
controller of the form (l2Tli that minimizes the cost function 

(Eoll. 

To obtain an explicit expression for J, we consider the 
closed loop system: 


dy = Ar/dt + Bwcl] 


where. 



WCL 


dwi 
dui ; 

dz;2 


dwcL dwcL^ = dt\ SwcL = ^^{Fwcl)'^ 


A — ^ BuCk 

ByC Ak + ByDuCK 

„ ^ \ By,, Bu O' 

ByDw^ ByDu + By^ By^ 

Finally, 

J = Tr {RQ) , (22) 

where Q is the unique symmetric positive definite solution of 
the Lyapunov equation 


AQ -f- QA^ + BSyjfy^B^ — 0 ; 


and 


Ri 0 
0 C^R2Ck 


For a more detailed derivation of these expressions, see 0. 

In the following subsection, we present an algorithm for 
designing a quantum controller of the form (l2Tl i which is a 
suboptimal solution to this problem. 


B. Design Algorithm 

We start by forming an Auxiliary Classical LQG Problem. 
Consider the plant and controller equations (fT9] l and (ISTT i 
and define du = du — dui. By temporarily ignoring the 
By^ dui and By^ du 2 noise terms, and treating dwi and dui as 
classical Wiener processes with intensity matrices Sy,^ and Sy.^ 
respectively, we obtain the auxiliary classical plant equations: 

dx = Ax dt + By du + By,^ dwi + By dui; 
dy = Cxdt -\- Dydu + Dyj^dwi + Dydvi] (23) 

and the auxiliary classical controller equations: 

dxK = Aj^xk dt + By dy; 
du = CkXk dt. (24) 


We also define an Auxiliary Cost Function which introduces 
an extra term to account for the fact that we have ignored 
the noise terms By^ dui and By^ dz ;2 that will appear in the 
quantum version of the controller: 


Jaux = lim sup —E 

tf—¥oo if 


x'^{t)Rix{t) + /3y(t)'^R2/3u(t) dt 


uo 


-I- lim sup —E 

tf—^oo if 


p/3„(f)^i?2/3u(f) dt 


(25) 


where E[.] denotes the classical expectation and p > 0. 
Equivalently, 


Jaux = lim sup —E 

tf—¥oo if 


x'^(t)Rix{t) + /3„(f)^i?2/3«(f)df 


L^o 


(26) 


where i ?2 = (1 + p)R 2 , p > 0. The Auxiliary LQG problem is 
to find a controller (l24l) that minimizes the cost function (l26l) 
for the plant (|2^ . 

Our approach to the coherent quantum LQG problem is 
as follows. The auxiliary LQG problem is first solved for 
a given p > 0 and the resulting auxiliary controller (l24l) 
is implemented as a quantum controller (EB by applying 
Theorem [3 or Theorem IH The cost function (l20l i is then 
evaluated using the expression (|2B . Linally, this process is 
repeated, optimizing the cost function (l20l) by using a line 
search over the parameter p to obtain our final suboptimal 
controller. 

We now detail one iteration of this design process. The 
auxiliary LQG problem is a standard classical LQG problem 
and is solved in the usual manner; see for example ll2^ . The 
solution is the auxiliary controller (l24l) with 


Ak = A-KC-ByF + KDyF; 
By = K; 

Ck = -F. 


Here, F and K are obtained as follows: 

F = R^^BlP; 

where P > 0 is the solution to the ARE 

Ri - PByk^^ bIp + A^P + PA = 0, 

















and 

K = {QC^ + Vi2)V^^-, 

where Q > 0 is the solution to the ARE 

0 = {A-Vi2V^^C)Q + Q{A-Vi2V^^CY 
- QC'^V^^CQ + El - ^12^2“ Vi^. 


Here, 


^Wi 

Bu 


dwi 


dwi 

T 

^Wi 

Bu 

Duji 

Du 


dui 


dui 



Du_ 


El Ei 2 
Vi^ ^2 


dt; 


El Ei2 

y^2 ^ 2 . 


Buui Bu 


'Su,, 0 ■ 


Bwi Bu 

Bwi Du 


. 0 S-K 


Duji Du 


Next, we obtain a fully quantum system of the form (|2TI) . 
based on the auxiliary controller (l24l) with , Sj,, Cif} 
obtained above. We first attempt to apply Theorem [3 If the 
conditions of the theorem are satisfied, the transfer function 
of the auxiliary controller is implemented as a system (ED, 
with only direct feedthrough quantum noises introduced by 
applying Theorem [3 That is, riy = n„. If the conditions 
of Theorem [3 are not satisfied, then the auxiliary controller 
is implemented by applying Theorem |3 which will result in 
> riu quantum noises. 

Finally, the cost function (l20l i is evaluated using the ex¬ 
pression (123) . For details on obtaining and By^ see 
Section IIV-AI 

Our algorithm is summarized as follows: 

1) For a given p > 0, form the Auxiliary Classical LQG 
Problem (l23 l. (l25l l. 

2) Solve to obtain the classical auxiliary controller (l24li . 

3) Implement this controller as a coherent quantum con¬ 
troller (I 2 D . 

4) Form the corresponding closed loop system, and evaluate 
the resulting cost function (l20l i. 

5) Repeat, optimizing over p > 0. 

We now give a heuristic motivation for our algorithm. In 
the standard (classical) LQG problem, the separation principle 
allows the optimal state estimator and optimal regulator to 
be designed independently and then combined to yield the 
optimal controller. In contrast to this, in the quantum version 
of the problem, the regulator gain Ck directly affects how 
strongly the quantum noises dui and dt ;2 impact the state 
estimator because By^ and By^ depend on Ck- 

Our method ignores the introduction of the additional noises 
dui and dt ;2 when designing our state estimator. In order to 
ensure that the effect of these noises is not too great, when 
designing the regulator we introduce the parameter p which 
puts an additional penalty on the size of the control signal. 
The final step of optimizing over p ensures the right balance: 
if p is too small the effect of the additional noises dui and 
dv 2 dominate the closed loop system response leading to poor 
performance whereas if p is too large, the feedback gain is 
unduly penalized also leading to poor performance. 


We justify our approach by observing that in practice, our 
algorithm is computationally tractable and examples show that 
the controllers so obtained yield good results. In particular, 
the example which follows demonstrates how a suboptimal 
coherent quantum controller can outperform a combination of 
heterodyne measurement and optimal (classical) measurement 
based feedback control. 

VI. Illustrative Example 
T o demonstrate the coherent quantum controller design 
process of Section |V] we consider a two mirror optical cavity 
driven by thermal noise of intensity as shown in Figure [3 
This example is a modification of an example considered in 
EH. Optically coupled to the second mirror is a controller to 
be designed via our algorithm. We compare the value of the 
cost function obtained with our controller to the no control 
case; i.e. when the second mirror is not connected to any other 
system and thus driven by a vacuum noise. We also compare 
our approach with a scheme involving heterodyne detection 
and an optimal classical LQG controller. Our design objective 
is to minimize the expected number of photons in the cavity: 
N = (a^a) = I -I- X 2 ^) — 0.5, where Xi and X 2 are 
the position and momentum operators for the cavity. It will 
be shown that for all > 0, it is possible to achieve better 
performance using a quantum controller designed using our 
method than with the optimal measurement based controller. 
This validates the utility of our method. 


Input 


Optical Cavity 


Controller 


Fig. 3. Quantum Plant: we wish to minimize the expected number of photons 
in the cavity. 


Our plant is of the form ( fT3 with 


A — — ^/2 x2 ; 


Bu 

= —\fi^l2x2', 

Byj-^ 

= s/^l2x2', 

c 

= s/^l2x2', 

Du 

= ^2 x 2; 

Dyy-^ 

= 02 x 2 ; 

Syj-i 

= (l + 2fc„)J2 

7 

= 0.2; 

Ki 

= 0.1; 

K2 

= 0.1. 


(27) 


It is sufficient to minimize the cost function ( l20b with Ri = 
12x2 and i ?2 = 02x2- Then N = jj - 0.5. 

Remark 6: We use R 2 = 10“®/2 x2 to design both the 
heterodyne classical LQG controller and the auxiliary LQG 
controller. However, the cost function ( l20b for the resulting 
controller is then evaluated using i ?2 = 02 x 2 - 






























A. No Control 


We first consider the case of no control as a reference. Here, 
the system input is driven by vacuum noise: du = dui. The 
resulting closed loop system is: 


dx = Ax df + \Bw-i 


Bu] 


dwi 
dui ’ 


where dwcL = [ ^ quantum Wiener process with 

covariance 




(1 + 2fc„)/2x2 0 

0 ^2x2 


i ?2 = 10 ® and hence {Ak, By, Ck} are obtained. The closed 
loop system is then as follows: 


da; 

da;^ 



dt + B dw; 


A 

B 


A BuCk 

ByC Ak + ByD^CK] ’ 
Bw 

ByD,^ 


The value of the cost function ( l20b can now be calculated 
using (I 22 I 1 with i ?2 = ^ 2 x 2 and substituting for Swcl- 


Equation (l22l i can then be used to hnd the value of the cost 
function (l20l) for a range of values for kn- 


B. Heterodyne measurement and classical LQG control 

Next we consider combining heterodyne measurement with 
a classical optimal LQG controller. Heterodyne measurement 
introduces an additional vacuum noise input. Similarly, the 
output of the classical controller will contain a vacuum noise 
component when applied to the input of the plant. This is 
accounted for by making the following substitutions: 

du = dti + dw3 ; 
dy = dj/ + dw2', 

into dull to obtain an augmented plant. Here du and dy are 
classical signals which represent the input and output of the 
augmented plant. The resulting equations for the augmented 
plant are as follows: 


da: = Axdt + Budu +B^dw] 
dy = Cx df + Du du + Z?™ dth. 


Here A, Bu, C, and Du are as before, and 

Bw = [— V^T2x2 02x2 — ^/^I2x2] 
Dw = [02x2 ^2x2 ^ 2 x 2 ] ; 


dw 


dwi 

dw2 ; 

dW3 


(1 + 2fc„)/2x2 0 0 

0 12x2 0 

0 0 ^2x2 


C. Quantum LQG control 

Finally we consider our proposed control scheme. First the 
auxiliary LQG problem is formed. The auxiliary plant is given 
by ([23]) with parameters as in ( |27] ). The cost function is given 
by (1251 ) with Ri = 12x2 and R 2 = 10“®. 

Then, optimizing J over p, we do the following: 

1) Solve the auxiliary LQG problem as detailed above to 
obtain {AK,By,CK}- 

2) Obtain a physically realizable quantum implementation 
of {AK,By,CK}- We do this by hrst attempting to 
apply Theorem 3 to obtain Bki- If the conditions of this 
theorem are not met we apply Theorem 2 to obtain Bki 
and Bk 2 - 

3) Evaluate the cost function (i20l i using (1221 ) with Ri = 
I 2 x 2 and i ?2 = 02 x 2 - 

D. Comparison of controller performance 

The relative performance of the no control case, the classical 
LQG case and our coherent control case are illustrated in 
Figure |4| and Figure [5] In the regime where both the ther¬ 
mal noise driving the system and the quantum noises are 
signihcant, the coherent quantum feedback controller offers 
the best performance of the schemes considered. If we leave 
the quantum regime with ^ 0, the relative performance 
benehts of the coherent quantum feedback controller over 
the measurement based feedback controller diminish as the 
thermal noise dominates the system and the quantum noises 
become insignificant by comparison. In the limit as 0, 

where the system is driven only by vacuum noise, our pro¬ 
posed controller offers no advantage over the no control case. 
This is consistent with the idea that the cavity cannot be driven 
below the vacuum state. 


As with the auxiliary LQG problem in our design algorithm, 
we treat dth as a standard classical Wiener process with 
intensity matrix S^. We now have a standard classical LQG 
problem. We wish to hnd a controller of the form: 

dxK = AkXk dt + By dy; 
du = CkXk dt. 

The estimator gain K and the regulator gain F are obtained 
in the usual manner. For this example R 2 = 02x2 but for 
computational reasons we assume R 2 takes a small value of 


VH. Conclusion 

The notion of physical realizability is fundamental to the 
coherent quantum feedback control problem where we wish to 
implement a given synthesized controller as a quantum system. 
By introducing additional quantum noises, it is always possible 
to make a given LTI system physically realizable. However, 
introducing additional quantum noises is undesirable in terms 
of the control system performance. In this paper, we have given 
an expression for the number of introduced quantum noises 
that are necessary to make a given LTI system physically 
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Fig. 4. Expected photon number N for no control, heterodyne measurement 
and classical LQG control, and coherent quantum LQG control. 


Normalized Photon Number: N / Nno control 



Fig. 5. Normalized photon number N/Njsic for heterodyne measurement 
and classical LQG control, and coherent quantum LQG control. 

realizable. Our result also gives a method for constructing the 
resulting fully quantum system. 

We also considered the case where a strictly proper transfer 
function is to be physically realized. We have given a result in 
terms of a Riccati equation for when it is possible to physically 
realize a specified transfer function by only introducing direct 
feedthrough vacuum noises and no additional quantum noises. 
We have also given conditions for when this Riccati equation 
has a suitable solution. 

Using these results we have developed an algorithm for 
obtaining a suboptimal solution to a coherent quantum LQG 
control problem. Our example demonstrates the utility of our 
results and shows that coherent quantum feedback control can 
offer performance benefits over measurement based feedback 
control. 
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